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Abstract 

Background: Long non-coding RNAs (IncRNAs) constitute a novel class of non-coding RNAs. LncRNAs regulate 
gene expression, thus having the possibility to modulate disease progression. In this study, we investigated the 
changes of IncRNAs expression in the heart after myocardial infarction (Ml). 

Results: Adult male C57/BL6 mice were subjected to coronary ligation or sham operation. In a derivation group of 
4 Ml and 4 sham-operated mice sacrificed 24 hours after surgery, microarray analysis showed that Ml was associated 
with up-regulation of 20 IncRNAs and down-regulation of 10 IncRNAs (fold-change >2). Among these, 2 IncRNAs, 
called myocardial infarction-associated transcript 1 (MIRT1) and 2 (MIRT2), showed robust up-regulation in the Ml 
group: 5-fold and 13-fold, respectively. Up-regulation of these 2 IncRNAs after Ml was confirmed by quantitative 
PCR in an independent validation group of 8 Ml and 8 sham-operated mice (9-fold and 16-fold for MIRT1 and 
MIRT2, P < 0.001). In a time-course analysis involving 21 additional Ml mice, the expression of both IncRNAs peaked 
24 hours after Ml and returned to baseline after 2 days. In situ hybridization revealed an up-regulation of MIRT1 
expression in the left ventricle of Ml mice. Expression of MIRT1 and MIRT2 correlated with the expression of multiple 
genes known to be involved in left ventricular remodeling. Mice with high level of expression of MIRT1 and MIRT2 had 
a preserved ejection fraction. 

Conclusion: Myocardial infarction induces important changes in the expression of IncRNAs in the heart. This study 
motivates further investigation of the role of IncRNAs in left ventricular remodeling. 

Keywords: Myocardial infarction, Left ventricular remodeling, Non-coding RNAs 



Background 

Cardiac diseases including stroke continue to be the 
main cause of death and disability in developed coun- 
tries [1,2]. Despite modern reperfusion strategies, a still 
significant proportion of patients develop left ventricular 
(LV) remodeling leading to heart failure after myocardial 
infarction (MI). Oxygen and nutrient deprivation to the 
heart induces severe damages, which can be of multiple 
types: necrosis or apoptosis of cardiac cells, cardiomyo- 
cyte hypertrophy, or fibrosis. Part of these damages can 
be induced by a de-regulation of gene expression. 

Since the initial sequencing of the human genome 
more than a decade ago [3,4], huge progress has been 
made in the understanding of its complexity. It appears 
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now that only a minor part of the human DNA encodes 
proteins, while the remaining is transcribed into non- 
protein coding RNAs [5-7]. Non-coding RNAs have 
been arbitrarily dichotomized as short non-coding RNAs 
(20-22 nucleotides -long, called microRNAs, miRNAs) 
and long non-coding RNAs (IncRNAs, >200 nucleo- 
tides). While miRNAs down-regulate gene expression 
mostly by destabilization of target messenger RNA [8], 
the regulation of gene expression by IncRNAs appears to 
be much more complex, involving both activation and 
repression of gene expression, and modulation of chro- 
matin architecture [9]. Since their discovery [10,11], 
IncRNAs have emerged as attracting biomarkers and 
therapeutic targets in the oncology field. However, our 
knowledge of the role of IncRNAs in cardiovascular dis- 
ease is only at its infancy. 

Only few studies reported associations between IncRNAs 
and the heart. Two landmark studies reported the 
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identification of 2 IncRNAs involved in cardiac devel- 
opment, Braveheart [12] and Fendrr [13]. In a genetic as- 
sociation study, the IncRNA MIAT (myocardial infarction- 
associated transcript) has been shown to be associated with 
the risk of MI [14]. Recent studies reported dysregulation 
of IncRNA expression in the failing heart [15,16]. However, 
the role of IncRNAs in the infarcted heart is still poorly 
characterized. In particular, whether IncRNAs may affect 
the course of LV remodeling post MI is unknown. 

The present study was designed to (1) determine the 
effect of MI on the expression of IncRNAs in the heart, 
and (2) identify IncRNAs potentially involved in LV re- 
modeling post MI. 

Methods 

Animal experiments 

This study was conducted in accordance with the regula- 
tions of the Animal Welfare Act of the National Institutes 
of Health Guide for the Care and Use of Laboratory Ani- 
mals (NIH Publication No.85-23, revised 1996). Protocols 
were approved by the Regional Veterinary Department 
('Direction Departementale de la Protection des Popula- 
tions'), agreements RAR1A03516811825 and 54-100. 

Mice were anesthetized by inhalation of isoflurane/ 
oxygen mixture (2.5%/1.5 v/v). When the mice were un- 
responsive to toe-pinch, they were intubated and venti- 
lated with a rodent respirator and were placed on a 
heating pad. A left thoracotomy of the third interrib 
space was performed to expose the heart. After pericar- 
dial incision, permanent occlusion of the anterior inter- 
ventricular artery was performed with a 7-0 Prolene 
suture. Having confirmed the presence of myocardial 
infarction by observation of ventricular blanching, ribs 
were closed with a 6.0 Vicryl suture, muscles were 
repositioned and the skin was sutured. The endo- 
tracheal tube was removed after spontaneous breath- 
ing. After surgery, mice were placed in an incubator at 
30°C for at least 30 min and then returned to their 
cages. 

For FDG-PET exam, mice received a pre-medication 
of 100 mg/kg of Acipimox in two intraperitoneal injec- 
tions. The first dose was injected 1 hour before the in- 
jection of FDG and the second dose simultaneously to 
FDG injection. Pre-medication with Acipimox allows a 
higher myocardial uptake of FDG and enhanced signal 
to noise and myocardial to blood activity ratios [17]. 
One hour before the exam, 37 MBq of FDG was injected 
in tail vein. Recording was performed during 40 min 
under continuous isoflurane anaesthesia using a dedicated 
small animal PET system (Inveon, Siemens, Knoxville, TN, 
USA). FDG uptake was determined on collapsed short-axis 
slices in each segment from the 17-segment LV division 
from the American Heart Association [18] using the QPS 
software [19]. Left ventricular (LV) end-diastolic volume 



(EDV), LV end-systolic volume (ESV) and LV ejection frac- 
tion (EF) were obtained from contiguous ECG-triggered 
short-axis slices using the QGS software. 

Mice were sacrificed with an overdose of isoflurane/ 
oxygen mixture. Blood was harvested by cardiac punc- 
ture. Heart was excised and immediately homogenized 
in Lysis Binding Buffer (mirVana isolation kit, Life tech- 
nologies) for extraction of total RNA using the mirVana 
isolation kit (Life technologies, Merelbeke, Belgium) ac- 
cording to manufacturer's instructions. On-column DNase 
I digestion (Qiagen, Venlo, The Netherlands) was per- 
formed to eliminate potential contamination with gen- 
omic DNA. Concentration and integrity of RNA were 
assessed using a Nanodrop spectrophotometer (Nanodrop 
products, Wilmington, USA) and a 2100 Bioanalyzer 
(Agilent technologies, Santa Clara, USA), respectively. 

Microarrays 

Agilent microarray platform was used with Low input 
Quick Amp labeling kit (Agilent) according to manufac- 
turer's instructions. Briefly, one-color spike mix was 
added to 200 ng of total RNA prior to amplification and 
labeling steps. Complementary RNA was purified with 
RNeasy Mini kit (Qiagen) and hybridized onto mouse 
SurePrint G3 microarrays (8x60K, Agilent). High-resolution 
microarray scanner (Agilent) and Feature Extraction soft- 
ware were used to scan the slide and extract raw microar- 
rays data. Microarray data are available in the NCBI 
Gene Expression Omnibus [http://www.ncbi.nlm.nih. 
gov/geo] under the accession number GSE46395. 

Pre-processing of raw data was performed using Limma 
[20] and VSN [21] packages rooted in the statistical com- 
puting environment R. Spots with a signal that was not 
significantly greater than the corresponding background - 
flag automatically established by Feature Extraction - 
were removed. Normalization between arrays was per- 
formed with Agilent spike-in probes. Principal component 
analysis was performed with R package ClassDiscovery 
[22]. After removing control probes and the probes that 
were detected in less than 4 samples per microarray slide, 
differentially expressed transcripts were determined using 
the i-test procedure within Significance Analysis of Micro- 
arrays version 3.09 which uses data permutations to esti- 
mate false discovery rate for multiple testing [23]. 100 
permutations were used in our analyses. Heatmaps were 
created using Cluster 3.0 and TreeView [24]. Functional an- 
notation and enrichment analysis of differentially expressed 
genes were performed with DAVID [25] (The Database for 
Annotation, Visualization and Integrated Discovery). 

LncRNA identification and construction of correlation 
networks 

Microarrays were re-annotated for IncRNAs on a probe 
level. Briefly, microarray probes without accession prefix 
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NM (i.e. protein-coding), according to the manufacturer, 
were aligned with IncRNAs from IncRNAdb [26], RefSeq 
[27] and Ensembl ncRNA [28] databases using BLAST + 
[29] rooted in the Perl environment. Only the probes that 
perfectly matched IncRNAs and did not match protein- 
coding RNAs (NM prefix transcripts from RefSeq) were 
considered as probes for IncRNAs. 

Spearman's rank correlation coefficient between se- 
lected IncRNAs and other remodeling-related tran- 
scripts from Gene database was obtained using the R 
package WGCNA [30]. This package computes the cor- 
relation coefficients and significance is determined by Stu- 
dent's test. Correlations with a P value < 0.05 were 
visualized under the form of a network using CytoS- 
cape (PMID: 21149340). 

Real-time quantitative PCR 

One ug of total RNA was reverse transcribed using 
Superscript II reverse transcriptase (Life technologies). 
Controls without reverse transcriptase were performed 
to ensure the absence of genomic DNA amplification 
during PCR. Real-time PCR was performed with IQ 
SYBR Green supermbc in a CFX96 apparatus (Bio-rad, 
Nazareth, Belgium). PCR primers were designed using 
the Beacon Designer software (Premier Biosoft, USA) 
(Additional file 1: Table SI). PCR conditions were as 
follows: 3 min at 95°C, 30 s at 95°C, and 1 min annealing- 
extension (40-fold). Optimal annealing-extension temperature 
was determined for each primer pair. The specificity of 
the PCR reaction was confirmed by melting curve ana- 
lysis. GAPDH was chosen as housekeeping gene for 
normalization. Expression levels were calculated by the 
relative quantification method (AACt) using the CFX 
Manager 2.1 software (Bio-Rad). Inter-run calibrator 
was used to normalize inter-run variations between 
separate real-time PCR runs. 

In situ hybridization 

Hearts were fixed in formalin during 24 hours, and em- 
bedded in paraffin. 5-um thick sections were performed. 
Expression of MIRT1 was assessed by in situ hybridization 
using the miRCURY LNA" microRNA ISH Optimization 
Kit (Exiqon, Vedbaek, Denmark) according to the manu- 
facturer's instructions. A scramble probe was used as 
negative control. Briefly, after deparaffinization with xy- 
lene and ethanol, sections were permeabilized with pro- 
teinase K (1 ug/mL). Then, sections were incubated with 
40nM double-DIG LNA™ MIRT1 probe (Exiqon) in 
hybridization solution (Sigma -Aldrich, Diegem, Belgium). 
Sections were washed and incubated with blocking solu- 
tion (Roche, Howald, Luxembourg), and then with sheep 
anti-DIG antibodies coupled to alkaline phosphatase 
(Roche). Revelation was performed with NBT-BCIP solu- 
tion (Roche) and the reaction was stopped with KTBT 



solution. Nuclei were stained with Nuclear Fast Red 
(Sigma- Aldrich) . 

In situ hybridization coupled to immunostaining 

To determine the cellular localization of MIRT1, in situ 
hybridization was performed as described above, with 
sheep anti-DIG antibodies coupled to fluorescein instead 
of alkaline phosphatase. Then, slides were subjected to im- 
munohistochemical staining with a rabbit polyclonal anti- 
body against sarcomeric alpha-actinin (Abeam, Cambridge, 
UK) to detect cardiomyocytes, a rabbit monoclonal anti- 
body against vimentin (Abeam) to detect fibroblasts, or a rat 
monoclonal antibody against CD45 (SantaCruz, Heidelberg, 
Germany) to detect leukocytes. Alexa Fluor" 635-coupled 
goat anti-rabbit antibody and Alexa Fluor 8 633-coupled 
goat anti-rat antibody were used as secondary antibodies 
(Invitrogen, Merelbeke, Belgium). Vectashield was used to 
reveal nuclei. Images were recorded on a confocal micro- 
scope (Zeiss Laser Scanning Microscope LSM 510 Carl 
Zeiss Microscopy, Oberkochen, Germany) using the 
LSM 510 META software (Carl Zeiss Microscopy, 
Oberkochen, Germany). 

Statistical analyses 

Results are presented as mean ± standard deviation (SD). 
Statistical analyses were performed with the SigmaPlot 
vll.O software. The Shapiro-Wilk normality test pre- 
ceded all analyses, t-test and Mann- Whitney test were 
used to compare two groups of continuous variables fol- 
lowing Gaussian and non-Gaussian distributions, respect- 
ively. Correlations between 2 variables were assessed using 
the Spearman test. Multiple group comparisons were per- 
formed using one-way analysis of variance and pairwise 
comparisons were performed using the Holm-Sidak 
method. All tests were two-tailed. A P value <0.05 was 
considered significant. 

Results 

Induction of Ml in mice - derivation group 

A derivation group of 8 mice, 4 subjected to MI through 
coronary ligation and 4 subjected to sham-operation, 
was used to study the effects of MI on the expression of 
IncRNAs in the heart using microarrays. 

Left ventricular function was assessed in all mice by 
FDG-PET 24 hours after surgery. As shown in Figure 1, 
this exam allowed the visualization of infarcted areas of 
the heart in MI mice, which were absent in sham- 
operated mice. Considering that infarcted areas display 
less than 50% FDG uptake (due to loss of metabolic ac- 
tivity of the cells, and to cell death), we calculated infarct 
size for each mice, which is the number of infarcted seg- 
ments related to the total number of 17 segments (see 
Methods for details). MI mice had in average 4 ± 2 in- 
farcted segments, which correspond to a mean infarct 
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Figure 1 FDG-PET exam in the derivation group. Cardiac function was assessed in sham-operated mice (n = 4) and Ml mice (n = 4) using FDG-PET, 
24 hours after surgery. A. Representative polar maps of a sham-operated mouse and a Ml mouse, in which the left ventricle is divided in 17 segments. In 
each segment, the indicated number corresponds to the percentage of FDG uptake. The segment was considered necrotic when the FDG uptake 
was inferior to 50%. Color scale shows the amount of FDG uptake, with darker color indicating lower FDG uptake (necrotic areas) and 
brighter color indicating high FDG uptake (viable areas). B. Vertical (up) and horizontal (bottom) long axis end-systolic PET images of a Ml 
mouse. White arrows indicate necrotic areas (absent or low FDG uptake). C. LV volumes, EF and infarct size as determined by FDG-PET image 
analysis. Infarct size was evaluated by the number of necrotic segments (<50% FDG uptake). 



size of 24% of the left ventricle. Sham-operated mice did 
not display any segment with less than 50% FDG uptake. 
LV end-diastolic and end-systolic volumes were signifi- 
cantly elevated, and the EF was reduced by 36% in aver- 
age in MI mice compared to sham-operated mice. 

Microarray experiments - derivation group 

Total RNA from cardiac tissue of the 8 mice of the deriv- 
ation group was extracted and used in microarray experi- 
ments. The Agilent SurePrint G3 Mouse Gene Expression 
microarray containing 55681 probes recognizing both 
coding and non-coding RNAs was used in these experi- 
ments. Microarray data were analyzed as described in 
Additional file 1: Figure SI. After pre-processing micro- 
array data with R packages Limma and vsn, and after re- 
moving control probes and the probes that were detected 



in less than 4 of 8 microarrays, 34049 probes were 
retained for further analysis. Principal component ana- 
lysis of these 34049 probes showed a clear discrimin- 
ation of sham-operated and MI mice, demonstrating 
that MI significantly affected the cardiac transcriptome 
(Figure 2A). Principal component 1 accounted for 42% 
of total variability and principal component 2 accounted 
for 22% of total variability. Then, genes differentially 
expressed between sham-operated and MI mice were 
identified using significance analysis of microarrays 
method. Genes with a fold-change in expression >2-fold 
and a q-value <5% were considered as differentially 
expressed. 704 genes (918 probes) satisfied these cri- 
teria (Figure 2B-C). 

Data mining using the DAVID database revealed that 
differentially expressed genes were highly involved in 
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Figure 2 Microarray analysis. Cardiac transcriptome of 4 sham-operated mice and 4 Ml mice (derivation group) was characterized 24 hours after 
surgery using microarrays. A. Principal component analysis showing the ability of gene expression data to discriminate Ml mice from sham-operated 
mice. B. M-A plot showing the distribution of the genes in the dataset. The vertical axis displays log2 transformed-fold change and the horizontal axis 
displays the average signal of each gene. C. Heatmap of differentially expressed genes. For B and C, red color indicates genes up-regulated 
in Ml mice and green color indicates genes down-regulated in Ml mice compared to sham mice. Black color indicates genes with comparable 
expression between Ml and sham mice. Significance threshold was 2-fold with a q-value <5%. 



inflammation-related pathways (such as Cytokine-cytokine 
receptor interaction pathway, Chemokine signaling path- 
way, and Toll-like receptor signaling pathway) (Additional 
file 1: Table S2). 

We developed an analytical pipeline to extract IncRNA 
data from microarrays (Figure 3A). From the 55681 
probes contained in the microarrays, we identified 26535 
(48%) probes corresponding to transcripts with NM 



prefix (i.e. mature messenger RNAs, mRNAs). From the 
remaining 29146 probes, we identified 2831 probes corre- 
sponding to IncRNAs represented in IncRNAdb, RefSeq 
or Ensembl ncRNA database. These 2831 probes corres- 
pond to 5% of all probes on the microarray (Figure 3B). 
From the 2831 probes, significance analysis of microarrays 
identified 20 IncRNAs significantly up-regulated in the MI 
group and 10 IncRNAs down-regulated in the MI group 
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Figure 3 Effect of Ml on IncRNAs expression in the heart. Microarrays performed with the 4 Ml and 4 sham-operated mice (Derivation group) 
were mined for IncRNAs data. A. Analytical pipeline used to identify microarray probes recognizing IncRNAs. B. Percentage of probes on the 
microarray corresponding to mRNA and IncRNA transcripts. C. Heat-map of IncRNAs differentially expressed between Ml and sham mice with a 
threshold fold-change of 2-fold and a q-value <5%. Red color indicates IncRNAs up-regulated in Ml mice and green color indicates IncRNAs 
down-regulated in Ml mice compared to sham mice. 



(Figure 3C and Table 1). Therefore, this discovery phase 
allowed the identification of 30 IncRNAs regulated 
after MI. 

Validation of the effect of Ml on IncRNAs 

We sought to confirm by quantitative PCR the differential 
expression of the top 10 IncRNAs identified in microarray 
experiments in the derivation group (5 up-regulated and 5 
down-regulated in MI mice compared to sham mice). 

First, we used the RNA samples of the 8 mice enrolled in 
microarray experiments (Derivation group). We observed 
that some IncRNAs, such as NR 002860, were expressed 
at a high level, whereas others, such as AK008868, were 
expressed at a low level (Figure 4A). Differential expression 
was confirmed for 8 of the top 10 IncRNAs. Differential ex- 
pression was not confirmed for AK008868, presumably 
due to its low level of expression and therefore low accur- 
acy of its quantification by PCR. For NR_001582, the dif- 
ferential expression was present albeit not statistically 
significant. Noteworthy, 2 IncRNAs displayed robust up- 
regulation in the MI group: NR_028427 (6-fold, P = 0.03) 
and ENSMUST00000100512 (12-fold, P = 0.03). This is 
consistent with microarray results in which these 2 IncRNAs 
were also the most highly up-regulated in MI mice (5-fold 
and 13-fold for NR_028427 and ENSMUST00000100512, 
respectively; Table 1). We named NR_028427 as Myo- 
cardial Infarction-Associated Transcript 1 (MIRT1) and 
ENSMUST00000100512 as MIRT2. 



Second, we measured the expression levels of the same 
top 10 IncRNAs in an independent validation group of 
16 mice (8 sham and 8 MI). The results obtained were 
very similar to that of the derivation group (Figure 4B). 
Indeed, differential expression was confirmed for 7 of 
the 10 IncRNAs. Statistical significance was lost for 
NR 038041, although the trend remained the same 
(slight down-regulation in the MI group). Most import- 
antly, MIRT1 and MIRT2 were, like in the derivation 
group, the most up-regulated IncRNAs in MI mice, with 
increases of 9-fold (P < 0.001) and 16-fold (P < 0.001), 
respectively. 

Third, 21 additional mice were enrolled in a time-course 
analysis. Mice were subjected to coronary ligation and 
sacrificed after 1, 3, 6, 16, 24, 48, and 72 hours (n = 3 per 
time-point). Hearts were harvested and expression of 
MIRT1 and MIRT2 was assessed by quantitative PCR 
(Figure 4C). Expression of MIRT1 and MIRT2 displayed 
progressive increases after induction of MI, reaching a 
maximum of 10-fold and 19-fold after 24 hours. 

Collectively, these results show that MI induces signifi- 
cant changes in IncRNAs expression in the heart. 

Localisation of MIRT1 in the heart 

Three additional mice were subjected to coronary ligation 
and 3 were sham-operated. The presence of MIRT1 in the 
heart 24 hours after surgery was revealed using in situ 
hybridization. As displayed in Figure 5 A, MIRT1 expression 
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Table 1 LncRNAs differentially expressed between Ml (n = 4) and sham-operated (n = 4) mice as determined by microar- 
rays using a fold-change >2 and a q-value <5% 



Agilent systematic name 


Agilent gene symbol 


Identified IncRNA 


Fold change 


q-value(%) 


ENSMUST00000100512 


Gm10872 


NR_045747 (Gm 10872) 


1354 


0.00 


NRJ328427 


583041 6P1 ORik 


NR_028427 (583041 6P1 ORik) 


524 


0.00 


NR_002860 


A130040M12Rik 


VL30 


2.71 


0.00 


chr2:1 581 69020-1 581 91 295_R 




N R_027887| N R_0 1 5463 (9430008G03 Ri k) 


2.64 


0.05 


NR_001582 


Speer5-ps1 


NR_027506|NR_001582 (Speer5-ps1) 


2.63 


0.00 


NR_026561 


Gm8884 


NR_026561 (Gm8884) 


2.56 


0.00 


NR_03071 1 


2210403K04Rik 


NR_030711 (Mir22hg) 


2.49 


0.26 


AK033703 


Hmga2-ps1 


NR_037996 (Hmga2-ps1) 


2.48 


0.00 


XRJ05817 


3300005D01 Rik 


N R_045080| N R_04508 1 1 N R_045079 (3300005 D0 1 Ri k) 


2.34 


0.00 


chr5:23 1 95569-2321 855 1 _F 




NR_038090 (AI506816) 


2.33 


0.16 


ENSMUST00000099574 


Gm10791 


NR_045889 (Gm10791) 


2.31 


0.00 


AK014779 


4833427F10Rik 


NR_045459 (4833427F 1 0Rik) 


2.30 


0.26 


chr5:1 00849748-1 00858532_R 




Adapt33 (NR_034038.1 Mouse 543041 6N02Rik) 


2.20 


0.00 


chrl 9:5771425-5848475_F 




NEAT1 


2.23 


3.43 


chr14:1 15443612-1 15445950_F 




NR_029382 (Mir17hg) 


2.16 


0.16 


NR_033578 


Gm 15645 


NR_033578 (Gm15645) 


2.13 


0.00 


NR_033596 


573041 6F02Rik 


NR_033596 (573041 6F02Rik) 


2.11 


0.05 


chr5:74487385-74497428_F 




NR_015531 (2700023E23Rik) 


2.09 


0.00 


NR_003368 


Pvtl 


NR_003368 (Pvtl) 


2.05 


0.05 


NR_028126 


6330407A03Rik 


NR_028126 (6330407A03Rik) 


2.01 


0.00 


ENSMUST00000161216 




ENSMUST000001 61216 (231 0075G1 7Rik) 


0.36 


1.82 


NR_015572 


1 81 0014B01 Rik 


NR_015572 (1810014B01 Rik) 


0.38 


1.82 


chrl 0:69669025-69686275_F 




NR_038041 (2310015B20Rik) 


0.39 


1.82 


AK008868 


221 0409 D07Rik 


NR_045360 (2210409D07Rik) 


0.42 


1.82 


chr4:1 1899775-1 1931875_F 




NR_040577 (1 7001 23M08Rik) 


0.46 


1.82 


NR_0 15464 


A330069E16Rik 


NR_0 15464 (A330069E16Rik) 


0.46 


1.82 


XR_1 04992 
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was up-regulated in the healthy part of the left ventricle of 
MI mice (remote area). In situ hybridization of MIRT1 
coupled with immunostaining of cardiomyocytes (sarco- 
meric alpha-actinin), fibroblasts (vimentin) and leukocytes 
(CD45) revealed that MIRT1 was expressed by fibroblasts 
(Figure 5B). 

Correlation between IncRNAs, infarct size and LV function 

We investigated the correlation between expression 
levels of the top 10 IncRNAs, infarct size and LV func- 
tion as assessed by FDG-PET in the 4 MI mice of the 
derivation group, 24 hours after coronary ligation. There 
was no statistically significant correlation at the 5% level, 



presumably due to the low number of mice (Table 2). 
Nevertheless, relatively strong negative associations were 
observed between expression levels of MIRT1, MIRT2, 
and infarct size (correlation coefficients of -0.80). On the 
other hand, positive associations could be noticed between 
MIRT1, MIRT2 and EF (correlation coefficients of 0.85 
and 0.80, respectively), suggesting that MIRT1 and MIRT2 
may beneficially affect LV remodeling post MI. 

Association between MIRT1, MIRT2, and LV remodeling 

To gain further insight into the potential role of MIRT1 
and MIRT2 into LV remodeling, we used the following 
approach. First, we identified a list of 47 genes known to 
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be involved in remodeling using the keywords ("myocar- 
dial infarction" AND "ventricular remodeling") AND 
("homo sapiens"[Organism] OR "mus musculus" [Organ- 
ism])" as a query to NCBI Gene database (="remodeling 
genes"). Among these 47 remodeling genes, 38 genes 



were detected on microarrays. Then, we used the micro- 
array data of the derivation group of 4 sham-operated 
and 4 MI mice to determine the correlation between the 
expression values of these 38 remodeling genes and the 
expression values of MIRT1 and MIRT2. This allowed 
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building networks of correlation between the IncRNAs 
and the 38 remodeling genes (Figure 6A). MIRT1 was 
significantly correlated with 18 remodeling genes, and 



MIRT2 was significantly correlated with 17 remodeling 
genes. Interestingly, while both IncRNAs were correlated 
with the same remodeling genes, Lif was only correlated 



Zangrando et al. BMC Genomics 2014, 15:460 
http://www.biomedcentral.com/1471-2164/15/460 



Page 10 of 14 



Table 2 Correlation between IncRNAs expression, infarct size and LV function as assessed by FDG-PET in 4 mice with 
Ml of the derivation group 





Infarct size 


EDV 


ESV 


EF 


ENSMUST00000161216 


-0.20 (0.92) 


0.38 (0.62) 


0.20 (0.80) 


-0.14 (0.86) 


NR_038041 


-0.80 (0.33) 


-0.07 (0.93) 


-0.57 (0.43) 


0.75 (0.25) 


NR_040577 


-0.40 (0.75) 


0.56 (0.44) 


0.01 (0.99) 


0.36 (0.64) 


NR_0 15572 


-0.20 (0.92) 


0.10 (0.90) 


-0.20 (0.80) 


0.27 (0.73) 


AK008868 


0(1) 


-0.95 (0.05) 


-0.69 (0.31) 


0.29 (0.71) 


NR_001582 


0.40 (0.75) 


-0.89 (0.11) 


-0.59 (0.41) 


0.17 (0.83) 


NR_002860 


-0.40 (0.75) 


0.76 (0.24) 


0.50 (0.50) 


-0.10 (0.90) 


NR_027887/NR_015463 


-0.20 (0.92) 


0.46 (0.54) 


0.19 (0.81) 


0.16 (0.84) 


NRJ)28427 (MIRT1) 


-0.80 (0.33) 


-0.48 (0.52) 


-0.75 (0.25) 


0.85 (0.15) 


ENSMUST00000100512 (MIRT2) 


-0.80 (0.33) 


0.60 (0.42) 


0.00 (1.00) 


0.80 (0.33) 



Correlation coefficients (P value) are indicated. 



with MIRT1 and Nos3 was only correlated with MIRT2 
(Figure 6A, right network). Figure 6B lists the 10 stron- 
gest correlations between IncRNAs and remodeling 
genes whose expression was regulated after MI. Corre- 
lations were highly robust, with coefficients close to 1 
for some gene pairs. Finally, we characterized the regu- 
lation of the expression of some remodeling genes in 
the group of 21 mice sacrificed at different time-points 
post MI (Figure 6C). All genes were up-regulated. We 
also measured Nppb, which encodes BNP, and we ob- 
served that, while most remodeling genes remained el- 
evated 3 days after MI, Nppb had returned to baseline 
levels after 48 hours. This kinetic was similar to that of 
MIRT1 and MIRT2. 

Together, these data support an association between 
the IncRNAs MIRT1 and MIRT2, and genes known to 
be involved in LV remodeling. 

Discussion 

In this study, we observed for the first time that MI in- 
duces a significant regulation of the expression of IncRNAs 
in the heart. Some of these IncRNAs were correlated with 
protein coding genes known to be involved in LV remodel- 
ing. These IncRNAs constitute novel candidates for future 
investigations of the therapeutic value of IncRNAs. 

Four groups of mice were used in this study. First, a der- 
ivation group of 8 mice was used to profile the expression 
of IncRNAs using microarrays. FDG-PET exams scheduled 
24 hours after induction of MI allowed to characterize in- 
farct size and LV function. In average, infarct covered one 
fourth of the left ventricle. EF was reduced by 36% 24 hours 
after induction of MI, which attests for a significant loss of 
LV function. End-diastolic and end-systolic volumes were 
increased, consistently with LV dilatation. Second, a valid- 
ation group of 8 independent mice was used to confirm 
microarray data. Third, a kinetic using 21 additional mice 
allowed characterizing the evolution of the expression of 



IncRNAs in the cardiac tissue after induction of MI. 
Fourth, 6 mice were used to study the localization of 
IncRNAs in the infarcted heart. This experimental design 
supports the robustness of our findings. 

A whole genome microarray was used for the discov- 
ery phase of our study. We observed that MI affected 
the expression of a significant number of genes (704), 
allowing a clear discrimination of MI mice from sham- 
operated mice. Many of these genes had a known link 
with the regulation of inflammation. This was expected 
considering that mice were sacrificed in the early inflam- 
matory phase that occurs in the first 24 hours post MI. 
Of note, this time-point was chosen to identify early 
triggers of LV remodeling which could be used to blunt 
or inhibit the development of LV remodeling at a very 
early stage after MI. 

An in-house analytical pipeline was developed to ex- 
tract IncRNAs data from microarrays. A similar ap- 
proach has already been used elsewhere [31]. We could 
identify 30 IncRNAs whose expression was regulated 
more than 2-fold and with a q-value <5% following MI. 
Of note, 12 of these 30 IncRNAs were also dysregulated 
in the heart of isoproterenol-treated mice [15]. We then 
focused on the 2 IncRNAs most differentially expressed 
between MI mice and sham-operated mice, MIRT1 and 
MIRT2. Up-regulation of these IncRNAs after MI was 
consistently observed in all groups of mice, and peaked 
after 24 hours. This up-regulation might be, at least in 
part, attributed to infiltration of inflammatory cells into 
the heart. This is supported by microarray data showing 
an up-regulation of the leukocyte marker CD45 in MI 
mice compared to sham mice (2.1 -fold) and of the 
monocyte/macro phage marker CD68 (1.6-fold). How- 
ever, CD45-positive leukocytes could not be detected by 
immunostaining in the remote area of the heart where 
MIRT1 expression is observed. This suggested that the 
increase of CD45 and CD68 measured by microarrays 
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(See figure on previous page.) 

Figure 6 Correlation between IncRNAs and remodeling genes. Microarray data from the derivation group of 4 sham-operated and 4 Ml mice 
were used in these analyses. A. Networks indicating the strength of the correlation between the IncRNAs MIRT1 and MIRT2 and 38 coding genes 
known to be involved in remodeling ("remodeling genes"). Remodeling genes differentially expressed between sham-operated and Ml mice are 
coloured, with darker colour indicating a strong differential expression. Red colour indicates a higher level of expression in Ml mice compared to 
sham-operated mice. Remodeling genes unaffected by Ml are in white circles. A q-value <5% was used as threshold for differential expression 
between sham and Ml mice (significance analysis of microarrays method). The thickness of the edges indicates the strength of the correlation between 
the IncRNAs and remodeling genes. Dotted lines indicate no correlation. A p-value <0.05 was used as significance threshold for correlation (Spearman's 
rank correlation and Student's test). B. List of the remodeling genes significantly (p < 0.05) correlated with IncRNAs and differentially expressed between 
sham-operated and Ml mice, r indicates correlation coefficient and p indicates p-value. C. Kinetic of the expression of remodeling genes after Ml. The 21 
mice sacrificed at different time-points after Ml were used in these analyses (n = 3 per time-point). P values obtained by ANOVA are indicated. 



come from inflammatory cells that are infiltrated in the in- 
farct lesion. In situ hybridization coupled to immunohisto- 
chemistry confirmed that MIRT1 is mainly expressed by 
fibroblasts within the remote area of the left ventricle. Fur- 
thermore, expression levels of MIRT1 and MIRT2 ap- 
peared to be negatively correlated with infarct size and 
positively correlated with EF. Although this is consistent 
with the known impact of infarct size on LV function, 
large infarcts inducing a deterioration of LV function with 
decreased EF, it also strengthens our assumption that ex- 
pression levels of MIRT1 and MIRT2 are not a mere con- 
sequence of inflammation. 

Inflammation is an important component of the re- 
modeling process. To address a potential link between 
the up-regulation of MIRT1 and MIRT2, and LV remod- 
eling, we used microarray data to determine the correla- 
tions between the expression values of genes known to 
be involved in LV remodeling, and the expression values 
of MIRT1 and MIRT2. Strong correlations were ob- 
served, suggesting that these 2 IncRNAs may function- 
ally regulate LV remodeling. However, this remains to be 
further explored and the contribution of these IncRNAs 
in the remodeling process further demonstrated. 

It would be tempting to investigate the expression of 
MIRT1 and MIRT2 IncRNAs in the human failing heart. 
However, there are no known homologs of these IncRNAs 
in human. 

In network analyses, both IncRNAs generally correlated 
with the same remodeling genes, except for Lif, which was 
only correlated with MIRT1 and Nos3, which was only 
correlated with MIRT2. This finding deserves further in- 
dependent validation. The possible interaction between 
MIRT1 and Lif on one hand, and between MIRT2 and 
Nos3 on the other hand, as well as its potential role in LV 
remodeling, needs to be further addressed. 

Nppb, the gene which encodes BNP, was correlated 
with MIRT1 and showed a 4-fold increase in expression 
after MI, peaking after 24 hours and returning to base- 
line levels after 48 hours. This kinetic was different from 
other tested remodeling genes, which all remained ele- 
vated until at least 72 hours after induction of MI. This 
observation merits to be validated but already points out 
a possible interaction between MIRT1 and BNP. 



Lgals3, which encodes the lectin galactoside-binding sol- 
uble 3, more commonly known as galectin-3, was signifi- 
cantly correlated with MIRT1 and MIRT2. This observation 
is relevant and in line with the role of galectin-3 in fibrosis 
and with its recently characterized value as biomarker of 
heart failure [32,33]. However, no evident relationship could 
be evidenced between circulating levels of galectin-3 and LV 
remodeling in survivors of acute MI [34]. 

Overall, our results suggest that IncRNAs may be in- 
volved in the regulation of several pathophysiological 
pathways leading to LV remodeling: inflammation (TNF), 
extracellular matrix turnover (MMP9), fibrosis (TGFB1 
and LGALS3), apoptosis (p53). 

Conclusion 

This hypothesis-generating study has led to the discov- 
ery of novel IncRNAs that may play functional roles in 
LV remodeling post MI. Further investigations are re- 
quired to demonstrate the therapeutic potential of these 
IncRNAs. 
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